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ABSTRACT 

We report the discovery of the resolved disk around HD 131835 and present the analysis and modeling of its 
thermal emission. HD 131835 is a ^15 Myr A2 star in the Scorpius-Centaurus OB association at a distance 
of 122.7;^}|g parsec. The extended disk has been detected to ^ 1.5" (200 AU) at 11.7 fim and 18.3 fim with 
T-ReCS on Gemini South. The disk is inclined at an angle of ^75° with the position angle of ^61°. The flux of 
HD 131835 system is 49.3 ±7.6 mJy and 84 ±45 mJy at 11.7 /im and 18.3 //m respectively. A model with three 
grain populations gives a satisfactory fit to both the spectral energy distribution and the images simultaneously. 

This best-fit model is composed of a hot continuous power-law disk and two rings. We characterized the grain 
temperature profile and found that the grains in all three populations are emitting at temperatures higher than 
blackbodies. In particular, the grains in the continuous disk are unusually warm; even when considering small 
graphite particles as the composition. 

Subject headings: circumstellar matter — infrared: stars — stars: individual (HD 131835) 


1. INTRODUCTION 

Planet formation and evolutionary history is imprinted on 
the nature and distribution of circumstellar debris. Circum¬ 
stellar debris disks are composed of dust and planetesimals 
that are in orbit around main sequence stars, analogous to the 
asteroid belt and Kuiper belt in our solar system. The materi¬ 
als in the disks are thought to be the sec ond generation where 
the du st is constantly being replenished (iBackman & Paresc^ 
Il993h fro m collisions between planetesimals and sublimation 
of comets (iHarper et al.lll984l) . Most of the debris disks are 
discovered through detecting their infrared excess and char¬ 
acterized by their spectral energy distributions (SEDs). How¬ 
ever, studying their SEDs alone provides only limited infor¬ 
mation. While the disk temperature and the total infrared (IR) 
flux can be well determined from the SEDs, the spatial in¬ 
formation of the systems is inaccessible. In addition, there is 
a possibility that the excess actually originates from an un¬ 
related foreground or background source contaminating the 
beam. 

In addition, with the SED information only, the location of 
the grains will be degenerate with the grain properties. For ex¬ 
ample, small grains further away from the star might be heated 
to the same temperature as large grains that are closer toward 
the center. Thus, a single-temperature debris disk does not 
necessarily indicate all the grains are located in a thin ring. 
A similar argument applies to the opposite scenario where 
small grains might be heated to a higher temperature com¬ 
pared to large grains at the same location. Therefore, for a 
multiple-temperature debris disk, the SED information alone 
is typically insufficient to distinguish between the cases where 
grains are located at the same spatial location (i.e., a ring) 
verses at multi ple spatially distinct locat ions (i.e., belts or ex¬ 
tended disks). iKennedv & WvattI (120141) argued that this de- 
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generacy is small for A-type stars, due to the truncation on 
the small end of the grain size distribution by radiation pres¬ 
sure. They proposed that most two-temperature disks around 
A-type stars probably arise from multiple spatial components. 
To confirm that the disk indeed has multiple spatial compo¬ 
nents, resolved images are required. 

If the disks can be spatially resolved, we will be able to 
characterize the spatial distribution of the grains and the ge¬ 
ometry of the disks. We can use distribution of grains to probe 
the underlying physics responsible for dust distribution, and 
constrain the locations and nature of the grain parent bodies. 
In addition, the morphology of the planetesimal belts might be 
closely related to the presence of unseen planets dynamically 
sculpting the disk. For example, as seen in mid-IR thermal 
emission, HR 4796A has an inclined ring with o n e lobe being 
r^5% b righter than the other (iTelesco et 5.1120001) . IWvatt et 
(119991) argued that the secular perturbations from a planet of 
mass > 0.1 Mj located close to the inner edge of the disk could 
introduce this brightness asymmetry. The planet would im¬ 
pose the forced eccentricity on the ring, causing one side of 
the ring being closer to the star and thus showing the "peri- 
center glow." Such incidence demonstrates the level of infor¬ 
mation provided from the direct imaging far surpasses from 
the SED alone. 

Here we have attempted to spatially resolve the disk around 
HD 131835, a system chary terized solely by its SED prior 
to this work. iRizzuto et al.l (120111) assign a 91% member¬ 
ship probability to the Upper Centaurus Lupus (UCL) mov¬ 
ing group (a subgroup of the Sco-Cen association) based 
on its Galactic location and velocity. HD 131835 is a 
young St ax (^15 Myr, based on the age of the UCL esti¬ 
mated bv lMamajek et al.l [20021: ^16 Myr, based on an anal- 
ysis of F-type pre -main-sequence members of t he group by 
iPecaut et al.ll2012l) with a spectral t ype of A2IV (iHoulOl 19821) 
at the distance of d = 122.7 pc (Ivan Leeuwenll2007l) . The 
infrared excess em ission of HD 131835 was first reported by 
iMoor et ^ (120061) through searching the Infrared Astronom¬ 
ical Satellite (IRAS) and Infrared Space Observatory (ISO) 
databases for the stars in the vicinity of t he Sun. The MIPS 
data for this source were first published in lChen et al.l (12012h 
which shows the MIPS 24 and 70 pm data for all of the 
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Ide Zeeuw et alJ (119991) ScoCen A-type stars and helps to put 
the excess emission for this star into context (e.g. it is one of 
only 4 UCL/Lower Centaurus Crux A-type stars with Ljr/L^ 
commensurate to beta Pic, > 10“^). IChen et aT] (120141) found 
that the SED for HD 131835 can be described using a two- 
temperature model. 

In contradiction to the traditional disk evolution scheme, 
CO gas in debris disks around A-type stars may be fairly com¬ 
mon. Typically, the primordial gas dissipates in a few million 
year time scale during the protoplanetary phase. By the time 
it transforms into a debris system, the CO gas level will be un¬ 
detectable. However, there are several ex ceptional cases. For 
exani ple, 49 Cet and HD 21997 are 40 (IZuckerman & Son^ 
l2012h and 30 (iTorres et al.l2008l) Mvr old A-type stars ho sting 
gas-rich disks (IZuckerman et al.lll995l : lMo6r et al.ll201 Ih . The 
gas is believed to be th e second generation, possiM y due to vi¬ 
olent comet collisions (IZuckerman & Song|[2012l) . In the case 
of HD 21997, there is an alternative explanation based on a re¬ 
cent ALMA observation which indi cates that the gas m ay be 
of primordial origin in this system (iKospal et al.l 120131) . An¬ 
other possible production mechanism involves constant resur¬ 
facing of the parent bodies and sublimation or photodesorp¬ 
tion of the CO ice. For the /3 Pic system, its CO distribution is 
particularly interesting because it is not the same as the dust, 
even more highly asymmetri c and may imply the presence of 
a second undetected planet (iDent et all 120141) . Not only has 
CO been detected and characterized for P Pic, 49 Cet, and 
HD 21997 but it has also been detected around 5 A-type stars 
in Upper Sco (Hughes 2014, private communication). 

Besides having the IR excess from the dusty debris, HD 
131835 hosts a de tectable amount of carbon monoxide gas 
(iMoor et all 1201 31 M oor et al. in prep). For HD 131835, 
iKastner et all (120101) first reported the nondetection of CO 
emission with the 30-m Institut de Radio Astrononi i e Mil- 
limetrique telescope. However, IZuckerman & Son^ (12012h 
argued that if the comet collision model is correct, then the 
H 2 /CO ratio is unconstrained and thus the upper mass limit 
of the CO ps fo r HD 131835 is 4.06 x 10-^ M^. Later, 
iMoor et ^ (120131) announced the discovery of detected sub¬ 
millimeter CO emission with the Atacama Pathfinder Exper¬ 
iment (APEX) radio telescope through a survey. Once the CO 
gas is well characterized, it can provide us valuable informa¬ 
tion on the disk environment and the dust-gas associations of 
such system with the coexistence of gas and debris at an age 
>10Myr. 

Here we report the discovery of the spatially resolved de¬ 
bris disk around HD 131835 at mid IR wavelengths, tracing 
the thermal emission from the grains. Currently, these images 
are the only resolved images of the system. In this paper, we 
present the analysis, modeling, and characterization of the de¬ 
bris disk around HD 131835. Our mid IR images are based on 
two epochs of observation as described in § [2l The data pro¬ 
cessing details, including image reduction, photometry, and 
PSF subtraction, are in § El In §ia we analyze the stellar 
properties and show that a three-component model is required 
to describe the system. In § [5l we characterize the grain tem¬ 
peratures in the system and investigate the possible grain com¬ 
positions. We also briefiy discuss how the current generation 
of imaging instruments and telescopes can improve our un¬ 
derstanding of debris disks in the near future. Finally, we 
summarize our modeling results in § [^ 

2. OBSERVATIONS 


TABLE 1 

Observations with Gemini S / T-ReCS 


Program ID 

PA 

UT Date 

Star 

Filter 

Integration 

(GS-20+) 

(deg) 




Time (sec) 

08A-Q-40 

55 

2008-05-08 

HD 131835 

Si-5 

318.5 





Qa 

1216.3 




HD 136422 

Si-5 

86.9 





Qa 

144.8 



2008-05-11 

HD 131835 

Si-5 

637.1 




HD 136422 

Si-5 

86.9 

09A-Q-19 

148 

2009-04-18 

HD 131835 

Qa 

3011.7 




HD 136422 

Qa 

724.0 



2009-04-19 

HD 131835 

Si-5 

3127.5 




HD 136422 

Si-5 

1013.5 



2009-05-23 

HD 131835 

Si-5 

1042.5 





Qa 

1505.8 




HD 136422 

Si-5 

434.4 





Qa 

434.4 



2009-07-12 

HD 131835 

Si-5 

1042.5 




HD 136422 

Si-5 

434.4 

Total 




HD 131835 

Si-5 

6168.2 





Qa 

5733.8 




HD 136422 

Si-5 

2056.0 





Qa 

1303.1 


We used the T-ReCS instrument on the Gemini South tele¬ 
scope to obtain mid-IR imaging of HD 131835 in two pro¬ 
grams, GS-2008A-Q-40 and GS-2009A-Q-19. As part of 
these programs, the star HD 136422 was also observed as a 
photometric and point-spread function (PSF) calibration stan¬ 
dard. The observations of the reference star were interleaved 
between the observations of the target. All the images are 
taken in the Si-5 (A^ = 11.66 yum, AA= 1.13 jim) and Qa (A^ = 
18.30 /im, AA= 1.51 /im) filters. In these bands, the T-ReCS 
pixel size is 0.08976" ± 0.00021" on the skyQ Due to the 
high sky background and instrumental thermal background, 
we used the standard mid-IR ABBA chop-and-nod observing 
technique with a 15" throw between chop-nod positions. 

Since the disk had not been resolved prior to the 2008 ob¬ 
servations, the first epoch data were obtained with an arbi¬ 
trary chop position angle. This turned out to be very close 
to the apparent position angle of the disk. To confirm that 
the feature was truly from the emission of the resolved disk 
and not from artifacts due to imperfect chopping and nod¬ 
ding, the second epoch of observation was obtained in 2009 
with the chop position angle roughly perpendicular to the pre¬ 
vious one. Observation conditions were generally very good, 
with diffraction rings visible in the calibration images. The 
sky was very transparent, with low water vapor (<lmm) mea¬ 
sured towards zenith for the most of the nights. The processed 
images of HD 136422 gave diffraction-limited resolution of 
^ 0.39" and ^ 0.54" at 11.7 jim and 18.3 jim respectively. 
HD 131835 (and HD 136422) were observed in two nights in 
2008 and four nights in 2009, with the total on source time 
of 6168 (2056) seconds at 11.7 /im and 5734 (1303) seconds 
at 18.3 /im. Details of the observations are summarized in 
TablelU 

3. DATA PROCESSING 
3.1. Image Reduction 

^ http://www.gemini.edu/sciops/instruments/t-recs/spectroscopy/detector 
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We reduced all the images following the standard mid- 
infrared data reduction procedures. We linearly combined the 
chop and nod frames within each image to remove the instru¬ 
mental and sky background. We then subtracted the average 
row background and corrected for the column offset due to 
small bias variations among different channels. The cores of 
the PSF standard star HD 136422 are fairly circular (as shown 
in Fig.[T]c and d) with the FWHM of ^ 0.39" and ^ 0.54" at 
11.7 fim and 18.3 fim respectively. The target images in both 
wavelengths show the extended and elongated emission com¬ 
pared to the PSFs. Before the PSF subtraction process, these 
mid-IR images already show prominent disk emission extend¬ 
ing beyond 1" along the semi-major axis in both bands. 

3.2. Photometry 

We performed aperture photometry on HD 131835 using 
HD 136422 as the photo metric reference. T he flux of HD 
136422 was taken from ICohen et ^ (119991) . We applied 
aperture correction to our photometry measurement on HD 
131835. Based on the observations of HD 136422, the PSF 
structure was fairly stable over the course of each night; the 
enclosed flux in an aperture radius of 1" usually varied only 
from 1 to 3 %. Thus, we adopted one relationship of aperture 
radius versus enclosed flux fraction for each night. Due to the 
nonuniform residual background, the enclosed flux fluctuated 
at large aperture radii. We characterized these fluctuations as 
the uncertainties of the radius-enclosed flux relationships. 

Next, we applied atmospheric extinction correction and 
color correction. The atmospheric extinction varied from 0.09 
to 0.19 magnitude per airmass among different nights in both 
bands. These values are consistent with th e extinction ob¬ 
served on Mauna Kea (iKrisciunas et al.lll987h . where the con¬ 
ditions are similar to Cerro Pachon. For color correction, 
since the 11.7 jim and 18.3 jim fell within the wavelength cov¬ 
erage of the IRS spectrum, we used its shape to compute the 
correction factors. Detailed information about the IRS data is 
provided in § [U The flnal flux of HD 131835 is 49.3 ± 7.6 
mJy for 11.7 jim and 84 ± 45 mJy for 18.3 jim. 

We have considered uncertainties from the process of mea¬ 
suring the target’s flux, the standard star’s flux, the aperture 
correction, and the extinction correction. However, by com¬ 
paring the photometry measurement from each image (33 im¬ 
ages at 11.7 yum and 14 images at 18.3 /im), we noticed the 
systematic errors dominated over the random errors. Thus, 
being conservative, we quoted the sample standard deviations 
of photometry measurements from the ensemble of images as 
the uncertainties of our mean fluxes. Our measured fluxes are 
consistent with the values from the IRS spectrum. Table[2]lists 
the flux measurements of HD 131835 from this work and from 
other published literature, and the SED is shown in Fig.[2l 

3.3. PSF Subtraction 

After linearly combined the ABBA frames within each tits 
image, we weighted combined all these individual image flies 
to produce one flnal high signal-to-noise ratio (SNR) image in 
each band. The uncertainty for each pixel in the flnal images 
was calculated by propagating the errors through the combi¬ 
nation process. The PSF was constructed by a similar process 
except for the weight used for stacking the target images. 

The target, PSF, and PSF-subtracted flnal images are shown 
in Fig. [T] The diagonal lines in the flgure represent the chop 
position angles. The chop position angle was rotated roughly 
90° from 2008 to 2009 due to the resolved disk structure as 


described in § [2l We scale the PSF to the target star flux, 
then subtract the stellar flux out from the target images. § 14.11 
describes the stellar flux characterization in detail. By pro¬ 
cessing the 2008 and 2009 data separately, we see that 2009 
observation alone confirms the detection of the structure and 
shows the position angle of the disk that is consistent with 
the observations made in 2008. Thus, we are confident that 
the elongated structure detected in 2008 is real, rather than an 
artifact arising due to imperfect chop-nod motion. 

In both bands, the central region of the disk is generally 
brighter than the outskirt. The detected SNR is greater than 6 
and 2 at 11.7 yum and 18.3 yum respectively. At 11.7 yum the 
extended structure is detected at 1 a out to ^ 1.5" and 0.7" 
in the major and minor axes, corresponding to ^ 180 AU and 
^ 85 AU from the star in projection. The 18.3 yum emission 
is resolved out to ^ 130 AU. Such extended disk structure in 
the mid-IR suggests the presence of warm dust far away from 
the star. 

4. MODELING AND ANALYSIS 

The infrared excess of HD 131835 has been measured 
in multiple wavelengths. Table [2] shows the summary of 
the measurements taken from the published li terature. HD 
13183 5 was detected in all four WISE bands (ICutri & et al.l 
l2012b . However, the WISE wl measurement represents a 
lower limit due to saturation, and w4 measurement is an up¬ 
per lim it due to source confu sion, as flagged in the WISE cat¬ 
alog bv ICutri & et al.l ( 20T^. The /RA5 p hotometry measure¬ 
ments are taken from Moor et all (l200d) . The S pitzer/MIFS 
24 yum and 70 pm measureni ents are taken from IChen et al.l 
dlOll and IChen et all dloT^. The object 145 6545-354138 
in the AKARI catalog (I Yamamura et al.l 120101) was identi¬ 
fied as HD 131835 bv iLiu et al.l (120141) . Among the four 
AKARI/FIS bands, onl y the source measured in 90 pm is 
confirmed and reliable (I Yamamura et al .1 120 1 Ol) . Finally, the 
longes t wavelength point o n the SED was measured in 870 
pm bv iNilsson et al.l (120101) . 

We apply color correction to goo d photometry mea sure- 
ments between 3 pm and 100 pm. iMoor et ^ (120061) had 
computed the color corrected fluxes for the IRAS points. For 
other points that fall within the wavelength coverage of the 
IRS spectrum, we use the shape of the IRS spectrum to de¬ 
termine the correction factors. The other photometry points 
are corrected based on the best estimated local SED shape. 
Color-corrected fluxes and references used to compute the 
correction factors are listed in Table [2l We do not consider 
correcting the 870 pm point since it is on the Rayleigh-Jeans 
tail so the correction factor should be close to unity. If a point 
is color-corrected, the corrected flux is used in the following 
modeling and analysis. Eig.[2]shows the SED (with the color- 
corrected fluxes). The SED shows a single hump infrared ex¬ 
cess with the flux peaking around 70 pm. The system has 
Lir/L* - 2 X 10 - 3 . 

The Spi t zer IR S data were obtained and analyzed by 
IChen et all dMTl . It is generally presumed that the MIPS 
fluxes will have a lower a bsolute calibration uncertainty. 
Therefore. IChen et ^ (120141) calibrated the IRS spectrum by 
tabulating the synthetic MIPS flux and comparing it to the 
measured MIPS flux. They presumed that the source did not 
vary with time and then scaled the IRS spectrum to match the 
MIPS 24 pm flux. More d etails of the spectra l extraction and 
calibration can be found in IChen et al.l (120141) . The calibrated 
IRS spectra do not show any obvious solid state emission or 
absorption features, suggesting HD 131835 is not a system 
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Fig. 1.— (a) and (b) are the images of HD 131835 from T-ReCS on Gemini South at 11.7 jum and at 18.3 jum. (c) and (d) are the corresponding images of the 
reference star HD 136422 scaled to the stellar flux of the target. Contours in (a)-(d) are spaced in logarithmic scales, (e) and (f) are the PSF subtracted images of 
HD 131835, and contours are spaced by the 1-a background noise level. In both wavelengths, the disk is resolved out to approximately 200 AU. The lines on the 
lower right-hand corners indicate the chop position angles for observations in 2008 and 2009 epochs. 


rich in small silicate grains and polycyclic aromatic hydrocar¬ 
bons. For the following analysis, we treat each point of the 
spectra as a single measurement that is directly comparable to 
photometry measurements. When including images in the fit, 
each pixel is considered as a data point and is treated equally 
as a point on the SED. For each image, only the central 68 
by 68 pixels 6" by 6") are considered in the fit. There 
are total of 9595 data points in which 10 points come from 
broadband photometry, 337 points come from the IRS spec¬ 
trum, and 9248 points from the images. 


4.1. Stellar Properties 

The measured fiux of the system includes the contribu¬ 
tion from the star and from the disk. The star has Tgff 
= 8770 K, log;^ = 4.0 and solar metallicity with A(V) = 
0.187 mag (IChen et al.ll2012h . We use the IDL Astro l ib rou - 
tine ccm_unred.pro, whic h is based onICardelli et al.l (Il989h . 
to apply reddening to the iKuruczl (Il993h stellar atmosphere 
model. Since the disk emission is prominent in longer wave¬ 
lengths, only the measurements with wavelength shorter than 
5.8 yum are used to fit for the stellar fiux to characterize the 
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TABLE2 

Summary of the Measured Fluxes oe HD 131835 


Wavelength (/xm) 

Source 

Flux (mJy) 

Color Corrected 

Flux (mJy) 

References 

0.43 

Hipparcos 

2571 ± 48 


(1) 

0.55 

Hipparcos 

2728 ± 26 


(1) 

1.24 

2MASS 

1448 ± 29 


(2) 

1.66 

2MASS 

965 ± 35 


(2) 

2.16 

2MASS 

652 ± 11 


(2) 

3.4 

WISE 

237.2 ± 5.8 


(3) 

4.6 

WISE 

166.8 ± 4.7 ^ 

165.3 ±4.7 

(3) (4) 

11.7 

Gemini S / T-ReCS 

49.7 ± 7.7 

49.3 ± 7.6 

(5) 

12 

WISE 

49.1 ±2.2^ 

49.4 ± 2.2 

(3) (4) 

18.3 

Gemini S / T-ReCS 

83 ±44 

84 ±45 

(5) 

22 

WISE 

160.5 ± 9.4 


(3) 

24 

Spitzer / MIPS 

153.1 ±3.1 

161.7 ±3.3 

(6) (7) (8) 

25 

IRAS 

186 ±34 

224 ± 40 

(9) 

60 

IRAS 

684 ± 62 

681 ±61 

(9) 

70 

Spitzer / MIPS 

659.2 ± 44.7 ^ 

710.0 ±48.1 

(7) (8) (10) 

90 

AKARI/EIS 

560 ± 39 

583 ±41 

(11) (12) 

870 

APEX/LABOCA 

8.5 ± 4.4 


(13) 

Spectrum 

5.2-37.9 

Spitzer / IRS 



(6) 


Reeerences. — ril lHc^g etiniMxil : 121 ICutri et allllOO^: Gl lCutri & r41 IWright et ( 5) this paper: 161 IChen et alllMTl 171 M IPS 

Instrument Handbook (Version 3.0, March 2011) by MIPS Instrument and MIPS Instrument Support Teams; 181 IRieke et al.ll2004l: 191 IMoor et al.ll2()0^: (10) 
IChen et al.ll20T^ : llll lYamamura et al.ll20TO : 1121 ILiu et al.ll20T^ : 1131 INilsson eU^IIoTO 


^ The quoted WISE data are based on aperture photometry. The error has included the uncertainty from the RMS scatter in the standard calibration stars. 

^ This is a lower limit due to saturation. 

^ This is an upper limit due to source confusion. 

^ The error is calculated by considering the source photon counting uncertainty, the detector repeatability uncertainty, and the absolute calibration uncertainty 



Fig. 2.— SED of the HD 131835 system. There are 15 well characterized 
photometry points and an IRS spectrum. The 11.7 /rm and 18.3 /^m T-ReCS 
photometry points are measured in this paper. Reference of other points are 
listed in Table|2] Note that the two triangle WISE points in 3.4 /rm and 22 /rm 
are the lower and upper limits. Well characterized photometry points between 
3 and 100 yum are color corrected. Only the points with wavelength < 5.8 /j,m 
are used to fit for the stellar atmosphere. The model is described in § 14.11 and 
the gray line shows the best fit. 


Stellar contribution. We compare the Hipparcos and 2MASS 
photometry to the synthetic fluxes of the model during the fit¬ 
ting process. The best fit is shown in Fig.[2l The stellar fluxes 
at 11.7 yum and 18.3 pm of this best-fit model were used to 
scal e the PSFs during the PSF subtraction process described 

in §[13 

We parameterize a scaling factor as ^ so that we 


can write 

nOO 

/ lydv. ( 1 ) 

Jo 

By integrating the unreddened stellar component of the best- 
fit SED model and solving the equation, we obtain ^ = 6.68 x 
10“^^. When estimating the luminosity, the largest uncertainty 
comes from the d istance measuremen t. Taking the distance to 
be 122.7±}|g pc (Ivan LeeuwerJ2007l) . we find/?* = 1.41 ±q {5 

Rq andL* = 10.5 ± 22 ^ 0 - 

4.2. Model 

Our goal is to And a disk model that can reproduce the ob¬ 
served SED and the images simultaneously. The SED sets 
strong constraints on the grain temperature distribution, where 
as the images inform the spatial distribution of the grains. We 
aim to recover the extended mid-IR emission, with the disk 
flux peaking close to the center of the system, while simulta¬ 
neously reproducing the broad IR excess hump on the SED. 
With the space and temperature distributions of the grains so 
constrained, we hope to further characterize the grain prop¬ 
erties such as its size and composition. In our model, we as¬ 
sume the disk to be optically thin since its Lir/L* ^2x 10“^. 
To search for a model that accounts for the extended emis¬ 
sion, we start with a simple two-dimensional continuous disk 
composed of only a single population of grains, then consider 
more complicated distributions. 

4.2.1. A Continuous Power-law Disk Model 

In contrast to perfect blackbodies, dust grains do not couple 
efficiently to radiation when the wavelength is much larger 
than the grain size. We consider a simple parameterization of 
the emissivity with a modified blackbody function, whereby 
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the emissivity varies with the frequency as a power-law with 
a positive index /3. The grain temperature Tg at a distance r 
away from the star is found by balancing the energy intake 
and the energy output: 

/V,* du = 471JBjy[Tg(r)] (2) 

here is the blackbody function. The grains are dis¬ 
tributed with the surface density of hq (r/lAf/)^, with hq be¬ 
ing the two dimensional number density of the grains at 1 AU. 
The density law applies to the region between the inner disk 
radius n and the outer radius it is zero elsewhere. The total 
disk flux /v,disk is 

r^o / X /3 / ^ \ r+i 

^ i (tTe) (tau) 

'i 

C is the scaling factor such that C = 2'K^nQa^{\Ps\])d~^. The 
characteristic grain radius is denoted as gq. 

There are seven free parameters total in this model. Five of 
which are needed to calculate the total disk flux: /3, T, C, n, 
and ro. Two additional parameters, inclination i and position 
angle 0, are needed to generate model images. We took sev¬ 
eral steps to initialize the starting parameters for the fit. First, 
we set P to zero as if the emission came from perfect black- 
bodies. Then, we set r, and to a range of a few pixels to 
several hundred AU. Afterwards, values for F and C were ob¬ 
tained by fitting the model to the SED data only. Finally, i and 
(j) were initialized to the best-fit inclination and position angle 
when fitting a ring to the images only. Once the parameters 
are initialized, we used a Levenberg-Marquardt least-squares 
fitter. The fitting process tends to drive to unbounded val¬ 
ues. The reason that tends to drift to unbounded values 
is due to the model’s under-prediction of the flux beyond 35 
/im. During the fitting process, after the model settles with 
the parameters that fit the majority of the data, the model will 
try to extend to an arbitrary large value so that it will con¬ 
tain slightly more cold dust. The contribution from this cold 
dust does make the SED fit better but the amount is negligible. 
Thus, we fix to 400 AU where beyond which the image fit 
does not get better. The best-fit parameters along with the to¬ 
tal chi-square (x^) and the reduced chi-square (xi) are listed 
in Table [3l The contributions to the total chi-square from the 
broad band photometry, IRS spectrum, and images are 3%, 
31%, and 66%. 

This one-component model does not produce a good fit. 
The modeled flux is too concentrated in the central region of 
the disk as the residual images show strongly negative values 
in center due to over subtraction. As a result, the best-fit i 
and (j) are not trustworthy. In addition, modeled SED under 
predicts the flux shorter than 16 /im and beyond 33 /im. The 
drive to better fit the short-wavelength part of the SED is in¬ 
consistent with the spatial location of flux in the images. This 
suggests that the power-law spatial distribution is inconsistent 
with the resolved outer disk flux and the contribution to the 
warm SED flux. Since this model clearly does not fit the data, 
we did not investigate further for estimating the uncertainties 
of the best-fit parameters. This inconsistency could be po¬ 
tentially resolved by considering a multi-component model. 
We therefore move on to a more sophisticated disk model by 
adding a second component. 

4.2.2. A Continuous Power-law Disk + A Ring 


This two grain populations model is composed of a contin¬ 
uous disk and a thin ring. The first component is as described 
in the section above. The second component assumes a single 
population of grains located in a narrow ring at a single ra¬ 
dius r away from the star. We adopt the same emissivity law 
in the form of for the ring component. However, this P 
parameter could have a different value than in the continuous 
component since the optical properties might be different for 
different grain populations. The flux from a ring model is: 

Fu,ringir) = CrBATgir)] (y^)^ (4) 

The grain temperature Tg is found by solving the Equation [2l 
and Cr is a scaling factor for the ring model. There are five 
parameters in the ring component: /3, C^, /, 0, and r. Lim¬ 
ited by the image resolution, we make the continuous disk 
and the ring sharing the same i and p. Here, we again fix 
Vo to 400 AU where beyond which the image fit does not get 
better and the SED fit does not show any significant improve¬ 
ment. We first use least squares fitting to And the best-fit 
model parameters. This model is a much better fit compar¬ 
ing to the previous single continuous disk model. The total 
decreases by more than 3000. The contributions to the 
total chi-square from the broad band photometry, IRS spec¬ 
trum, and images are 3%, 14%, and 83%. The image fit 
improves significantly. The residual images do not suffer 
from the severe over subtraction in the central regions any¬ 
more. Therefore, it is worth exploring the uncertainties of 
the fitting parameters for this model. However, the covari¬ 
ance matrix from the least squares fitting does not always 
provide reasonable error estimations, especially for parame¬ 
ters that are degenerate. Degeneracies are prominent between 
P, F, C, and n for the continuous disk model and between 
P, Cr, and r for the ring model. For example, effects of tem¬ 
perature depression from increasing r can be compensated by 
having a larger value of p. To better quantify the uncertain- 
ti es on the parameters , we us e the ensemble MCMC method 
of iGoodman & We^ (I201QI), as implemented in the emcee 
python package jForeman-Mackev et ^l2013h . 

The emcee package is based on the M arkov chain Monte- 
Carlo (MCMC) method developed by iGoodman & Wearel 
(I201QI) where they utilize multiple "walkers" to propagate 
multiple chains simultaneously. We set the initial positions of 
the walkers by drawing random numbers based on Gaussian 
distributions with means equal to the best-fit parameters and 
the variances equal to 10% of the means. This MCMC algo¬ 
rithm adjusts the candidate step proposal distribution based on 
the walkers’ current positions in the parameter space. We use 
100 walkers and run for few thousand steps after the burn-in 
period seems to be over. All of the marginalized probability 
density functions (PDFs) look fairly symmetric. We calcu¬ 
late the uncertainty on each parameter by measuring the la 
interval on the marginalized PDFs, with the upper and lower 
bounds measured separately. The parameter values and their 
uncertainties are listed in Tabled 

The best fit shows that the continuous component contains 
hotter grains and is mainly contributing to the 11.7 /im emis¬ 
sion while the ring component contains cooler grains and 
dominates the 18.3 /im emission. The hot continuous disk 
extends from a few tens of AU to hundreds of AU and con¬ 
tains grains as hot as 300 K. The cool ring locates around 
60 AU, spatially overlapping with the continuous component. 
By comparing to the single-disk model, this two grain popula- 
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TABLES 

One-, Two-, and Three-population Model Fits 


Parameter 

A Continuous Disk 

A Continuous Disk + A Ring 

A Continuous Disk + Two Rings 

Unit 

xl 

9920 

1.03 

6255 

0.65 

5958 

0.62 


Continuous Component 

f3 

0.76 

1.53!«“ 

1.64 ±0.02 


C 

3.7 X 10-31 

(5 ± 1) X 10-42 

(9±2)x lO-^'* 

[AU-i] 

F 

1.0 

0.48 ±0.05 

0 53+1111^ 


i 

49 

73.8±1.0 

74 5+11-9 

[deg] 


74 

62.0+}^ 

61 2^1-^ 

[deg] 

G 

0.24 

37±3 

35±3 

[AU] 

ro 

400 

400 

310'^^o 

.310_2o 

[AU] 

Ring Component I 

/3 


0 27^113 

0.59 ±0.02 


Cr 


(7.4 ± 1.5) X 10“‘5 

X 10-16 


r 


61^8 

105 ±5 

[AU] 

Ring Component II 

/3 

Cr 

r 



0.32 ±0.06 
(4±l)x 10-14 

220 ±40 

[AU] 


Note. — The chi-square (x^) and the reduced chi-square (x^) are the minimum values estimated by least squares fitting. The one-component parameters 
are from the least squares fitting; uncertainties are not estimated since the model does not describe this system well at all. The two- and three-component model 
parameters and uncertainties are quoted from the marginal distributions of MCMC results. For the two-component model, the parameters (from top to bottom) 
correspond to the quoted least squares are (/3, C, F, i, 0, ri, ro, /3, Cr, r) = (1.53, 5.0 x 10~^^, 0.47, 74.0, 61.8, 37.5, 400, 0.27, 7.2 x 10“^^, 62) For the three- 
component model, the parameters correspond (from top to bottom) to the quoted least squares are (/3, C, F, i, cj), ri, ro, (3, Cr, r, (3, Cr, r) = (1.64, Fix 10“^^, 
0.48, 74.3, 61.1, 36, 316, 0.59, 6.2 x lO-^^ 105, 0.33, 4.3 x lO-^^, 232). 


tions model significantly improved the fit to the images and to 
the short wavelength part of the SED. However, this model is 
still under predicting the fiux longer than 35 jim on the SED. 
The main reason for causing this underprediction is the out¬ 
numbered IRS spectral data points. The fit is driven by closely 
matching the spectrum in order to effectively lower the As 
a result, the model completely misses the photometry points 
at wavelengths longer than 35 /im. As a test, we considered 
excluding long wavelength IRS points. We find that in or¬ 
der to produce a reasonably good fit with this two-component 
model, we have to cut out IRS points longer than ^15 /im. 
Since we do not have a physical explanation to place such cut, 
we include all IRS points in our analysis. To account for the 
long wavelength emissions, we add another ring component 
in the following section. 

4.2.3. A Continuous Power-law Disk + Two Rings 

This three grain populations model is composed of a con¬ 
tinuous disk and two thin rings. The structures of these model 
components are described in the previous two sections. Al¬ 
though the inclination and the disk position angle are free pa¬ 
rameters in the model, we constrain the three components to 
share the same values. As before, we use the least squares 
fitting to find the best-fit model parameters and MCMC to 
quantify their uncertainties. The model parameters are listed 
in Table [3l The contributions to the total chi-square from the 
broad band photometry, IRS spectrum, and images are 2%, 
12%, and 86%. Eig.[3]and Eig.[4]show the best-fit SED model 
and the image models respectively. 

This best-fit model is composed of a hot continuous power- 



Wavelength (/xm) 


Fig. 3.— SED of the debris disk around HD 131835 with the best-fit three- 
population model (solid line) and its components (dashed lines). The residu¬ 
als are plotted in the bottom panel. This model is composed of a hot continu¬ 
ous power-law disk, a warm ring, and a cold ring. The grains are assumed to 
emit like modified blackbodies such that the emissivity is proportional io . 
The best-fit parameters are listed in Tabled Given this grain emissivity law, 
this three-component model is the simplest model that provides a reasonable 
fit to the SED while fitting the images simultaneously. 


law disk extended from 35 ± 3 AU to 310;^2o ^ 
perature ranging from 330 K to 150 K, a warm ring with 
temperature of 97 K located at 105 ± 5 AU and the cold ring 
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(a) 11.7 jim 
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0.075 

0.060 
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0.030 

0.015 

0.000 
mJy 
0.105 

0.090 

0.075 

0.060 

0.045 
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t>0 
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a 



1.5 1.0 0.5 0.0 - 0 . 5 - 1 . 0 - 1.5 


a 



1.5 1.0 0.5 0.0 - 0 . 5 - 1 . 0 - 1.5 


Aa(") 


Aa(") 


Fig. 4.— Images showing the best fit of the three-component model (a continuous power-law disk + two rings). Top panel: the PSF subtracted images, 
smoothed to suppress high-spatial-frequency noise for the displaying purpose. We smoothed the images with a gaussian with its FWHM (yellow circle) equals 
to 0.1". Contours are spaced by the l-cr background noise level derived from the smoothed variance map. Middle panel: the best-fit model convolved with the 
stellar PSF. Bottom panel: the smoothed residual flux divided by the uncertainty derived from the smoothed variance map. Images are smoothed with gaussians 
with FWHM equal to their PSF sizes only for the displaying purpose. Comparing large spatial structures in the residual images, this model seems to slightly 
under predict the central flux. Nonetheless, the variations in the resolution-element scale are within 1 (j, suggesting the model fits the images reasonably well. 
The best-fit parameters are listed in Table [3] 


with temperature of 52 K located at 220 ± 40 AU. The middle 
panel of Fig. [4| shows the PSF convolved model disk in the 
two imaging wavelengths. Approximately 90% of the flux 
in the 11.7 jim model comes from the hot continuous com¬ 
ponent while 10% comes from the warm ring. At 18.3 /im, 
about 2/3 of the flux comes from the warm ring and 1/3 of 
the flux comes from the hot disk. This model provides a rea¬ 
sonable lit to the images, despite there are some small patches 
of over-subtraction and under-subtraction as seen from the im¬ 
age residuals (the bottom panel in Fig.®. We characterize the 
variations in the resolution-element scales by dividing the sum 
of the smoothed residual flux values in a resolution-element¬ 


sized area by the sum of the corresponding smoothed noise. 
The variations in the resolution-element scales are within 1 
cr, suggesting the model fits the images reasonably well. The 
cold ring component is too faint to affect either imaging chan¬ 
nel. However, this cold ring is essential for contributing to the 
long wavelength part of the SED and is a key reason why we 
can constrain r^, unlike the previous two models. We found 
that this three-component model is be able to provide a sensi¬ 
ble fit to the SED and the images simultaneously. 

4.2.4. Other Disk Models 
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We also tried using a broken-power-law disk model. This 
model is composed of one continuous-disk described by two 
power-law grain distributions: one between the inner radius 
and the intermediate radius and the other between the inter¬ 
mediate radius and the outer radius. The entire disk shares 
a single value of emissivity power index /3, in which we as¬ 
sume the grain composition is the same across the disk. The 
best-fit result improves slightly compared to the one from a 
single continuous disk model described by only one power- 
law grain distribution (§ 14.2.11) . However, this best-fit broken- 
power-law model suffers from the similar inconsistency as the 
single-power-law model where the model fiux is too concen¬ 
trated in the center of the disk in both bands. The reason lies 
under the assumption of a single grain population. Since the 
temperature of the entire disk is governed by a single emissiv¬ 
ity power index, grains further away will have lower temper¬ 
atures. Thus, although a steeply rising density power-law can 
make the model disk fiux more extended, it will also drive the 
corresponding SED model too bright at longer wavelengths 
since the flux is originated from low-temperature grains. In 
order to have warm dust grains be responsible for the ex¬ 
tended emission instead, we must introduce a sec ond emis¬ 
sivity power index to our model (as used in § 14.2.21) . 

Another model we tried is composed of two continuous 
disks. We found that the two components can share the same 
values for T, r/, r^, /, and (j) while still producing a good 
fit compared to having two independent sets of parameters. 
However, this model is very sensitive to the p values. A small 
adjustment in /3 will introduce a significant deviation in the 
resulting model. Thus, we must keep the two /3 parameters 
independent in order to maintain a good fit. This result again 
indicates there are more than one population of grains in the 
system. We notice that one of the continuous components 
can be simplified into a ring without changing the best-fit re¬ 
sult significantly (xj, changes from 0.65 to 0.64). Therefore, 
we performed the detaile d analy sis of a continuous power-law 
disk plus a ring model (§ 14.2.21) instead of the two-continuous 
component model. 

5. DISCUSSION 

Our models support the inference of multi-spatial 
components from a m ulti-temperature SED. Recently, 
[Kennedy & WyattI (l2014l) argued that if the SED of a debris 
disk around an A-type star shows multiple temperature com¬ 
ponents, it is an indication that the system hosts multiple pop¬ 
ulations of grains in different locations. Our SED modeling of 
the disk around HD 131835 indicates there are multiple tem¬ 
peratures. Erom modeling with the resolved images in mid- 
infrared, we confirm the system indeed have grains at multiple 
spatial locations: a hot continuous component extends from 
35 ± 3 to 310:^20 ^ warm ring located at 105 ± 5 AU, and 

a cold ring located at 220 ± 40 AU. Our model indicates the 
two separated narrow rings are embedded in an extended disk 
component. Although not all the model components are com¬ 
pletely spatially separated, we are confident that the dust is 
not concentrated in a singl e belt. Our modeling resu lt agrees 
with the argument made bv iKennedy & Wv^ (l2014l) . adding 
to the small poll of high resolution observations confirmed 
cases. 

Grains with effective temperatures hotter than blackbodies 
are responsible for the observed disk emission, since perfect 
blackbody grains at these spatial locations would not have the 
appropriate color temperatures. In our models, the grains are 


assumed to emit according to a modified blackbody function 
where the emissivity is proportional to . In the case of 
the three-population model, the grains in the [3 values for the 
continuous component, the warm ring, and the cold ring are 
1.64 ± 0.02, 0.59 ± 0.02, and 0.32 ± 0.06 respectively. /3 = 0 
corresponds to blackbody grains. The positive P values indi¬ 
cate the grains are small since they are inefficient in absorbing 
and emitting at long wavelengths, and /3 = 1 corresponds to 
small grains in the limit lira A, where a is the grain size 
and A is the observed wavelength. Eor /3 > 1, real materi¬ 
als with complicated emissivities are required to explain the 
observations. Observational evidence points that it is com¬ 
mon for debris d isks around A-star s to have p values approxi¬ 
mately within 2 (iBooth et al.ll2013h . and our results fall within 
this region. 

We searched for the potential grain compositions by first 
examining the equilibrium grain temperatures with our modi¬ 
fied blackbody models and then compare them to the equilib¬ 
rium temperatures of specific grain compositions. Using the 
model parameter r and p, we plotted the corresponding equi¬ 
librium temperatures at different stellar locations in Eig. [5l 
Our goal is to use these temperature curves to identify the pos¬ 
sible grain compositions. We consider some common compo¬ 
sitions such as graphite, astro-silicate, and ic e. We use the 
optica l co nstants of graph i te and silicate from iDraine & Leel 
(119841) and lLaor & Drainel (119931) . The optical property of ice 
is a f unction of its t emperature due to varying crystal struc¬ 
tures. IWarrenI (119841) provided the index of refraction of pure 
ice at -1°C, -5°C, -20°C, and -60°C. Although the opti¬ 
cal properties are different at different crystal temperatures, 
the variation is quite small. Thus, to first order, we adopt 
the optical constants of ice assuming the ice preserves the 
similar crystal structure as when it is at T = -60°C. We cal- 
culated their equilibrium t emperatures using the Mie theory 
(iBohren & Huffma^ll983l) with the assumption that grains 
are spherical. The grain size contours are overplotted on the 
temperature-stellocentric distance map on Eig.[5l 

Grains composed solely of graphite, astro-silicate, or ice 
in the equilibrium temperature can not be used to explain all 
the observed disk properties. Among the three compositions, 
graphite grains have the highest temperature, and ice grains 
are the coolest. We have also considered using the other types 
of silicates such as olivine and pyroxene. However, their op¬ 
tical properties are similar to the astro-silicate, and thus, the 
grain temperatures are on the same order of magnitude as the 
astro-silicate grains. The warm and cold rings could be made 
up by graphite and silicate grains closer to micron size. On the 
other hand, the hot continuous disk component would require 
the emitters to be really small. Since small grains become 
inefficient emitters at long wavelengths, these grains will sus¬ 
tain much higher temperatures than a blackbody. However, 
even though nanometer-sized graphite grains are significantly 
hotter than blackbodies, they still could not match the temper¬ 
ature of the hot continuous disk. In addition, grains this small 
are not likely to be present in the system as discussed below. 

We have qualitatively considered having porous grains and 
using a broken emissivity power-law in attempt to address the 
abnormally warm dust. However, neither considerations are 
in our favor. Poro us dust grains have lower tem peratures than 
compact spheres (iKirchschlager & Woljl2Q13l) . Thus, intro¬ 
ducing porosity will make the discrepancy worse for finding 
a physical grain composition for the observed dust temper¬ 
atures. Similar qualitative result applies for using a broken 
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Fig. 5.— Equilibrium temperatures of the disk model components in com¬ 
parison to various grain types. The temperatures and radii for the three- 
component model are shown here. The red curve represents the hot contin¬ 
uous disk component. The green and blue dots represent the warm and cold 
ring components respectively. The black line is the temperature curve assum¬ 
ing grains are perfect blackbodies. Note that all three populations of grains 
have equilibrium temperatures higher than blackbodies. This grain property 
is commonly observed since small grains are inefficient in emitting at long 
wavelengths. The grain size (grain radii in /j,m) contours for graphite, sili¬ 
cate, and ice are overplotted here for comparison. Although the optical prop¬ 
erty of ice changes with temperature, we do not take it into account since the 
variation is quite small; here we use the T = 213 K ice crystal structure. Al¬ 
though small graphite grains, as expected, are significantly hotter than other 
grains, it still can not match the temperature gradient for the hot continuous 
disk component. Thus, grains composed solely of graphite, astro-silicate, or 
ice in the equilibrium temperature can not be used to explain all the observed 
disk properties. 


emissivity power-law case. A broken emissivity power-law 
assumes that the emissivity is 1 when v > Uc and 
when zz < z/c, where Vc is a free parameter indicating the criti¬ 
cal frequency at which the emissivity function changes. Com¬ 
pared to the smooth emissivity law we used, the grains with 
the broken emissivity power-law are less efficient in absorb¬ 
ing the stellar light. Therefore, grains with a broken emis¬ 
sivity power-law will sustain lower temperatures. In addition 
to making the grains cooler, a broken emissivity power-law 
introduces an additional model parameter, z/^, that makes the 
model more complex than a smooth emissivity power-law. 

Although the high disk-temperature draws the connection 
to small grains, no features of small grains are detected in 
the system. Nanometer-sized grains are in the molecular re¬ 
gion. For HD 131835, we do not observe any solid-state fea¬ 
tures in the 10 and 20 micron regions from the IRS spectrum, 
which indicates that the grains are probably la rger and more 
likely to be in the micrometer size region (e.g.. lDraine & Le3 
11984]) . Stochastically heating small grains such as polycyclic 
aromatic hydrocarbons (PAHs) could be one way to make 
grains have much higher t emperatures than in the equilibrium 
states (iDraine & Lill2001h . However, the IRS spectrum does 
not show the spectral signatures of PAHs either. 

Furthermore, small grains are unlikely to present in the sys¬ 
tem when considering the radiation pressure from the star. 
The blowout grain size is calculated by balancing the radia¬ 
tion pressur e with the gravity . Taking the mass of the star to 
be 1.9 Mq (IChen et al.fl2012b and using the grain density of 
p = 3.5 g/ crd , we get the blowout grain radius for HD 131835 


to be: 


^blowout — 0.91 pm 


10.5 L©; 


1^ 1-9 Me ^ 1^ 3.5 g/crn^ ^ 


(5) 

Grains smaller than this blowout size are unlikely to survive 
in the system. Although small grains can be trapped in reso¬ 
nance due to gas drag while migrating, the effect is unlikely to 
be significant with the gas level present in this system. Since 
small grains are unlikely to be responsible for the hot con¬ 
tinuous disk emission, the nature of these abnormally warm 
grains is not completely understood. 

Identifying unique disk models around A stars can be chal¬ 
lenging. For example , recently, with Her schel observations in 
70, 100, and 160 pm. lBooth et al.l (12013h shows that the disks 
around A stars have various morphologies, ranging from sys¬ 
tems that can be fit with just a narrow ring to the ones that 
require wider or multiple rings. Although direct images can 
provide constraints on grain size distribution and dust prop¬ 
erties, sometimes finding a proper model can be quite diffi¬ 
cult. For example, a detailed modeling work on [3 Leo de¬ 
bris disc with multi-wavelength observations shows the de- 
generacy between on e, two, and three components models 
(IChurcher et aLlUOlU) . For HD 131835, the three-component 
model gives a reasonable fit but more data are needed to set 
better constraints. We would like to better characterize the 
grain properties, understand the distribution of grains, and 
constrain the location of planetesimal belts in much greater 
detail. 

Fortunately, with the new generations of high resolution, 
high sensitivity, and high contrast instruments, detailed disk 
characterization is foreseeable. HD 131835 is located in the 
Southern sky, making it a favorable target for GPI (Gem¬ 
ini Planet Imager), SPHERE (Spectro-Polarimetric High- 
contrast Exoplanet Research), and ALMA (Atacama Large 
Millimeter/sub-millimeter Array). GPI and SPHERE can po¬ 
tentially image the disk in scattered light in the near infrared, 
providing the currently uncharacterized scattering properties 
of the dust grains around HD 131835, in both total inten¬ 
sity and linearly polarized light. Designed specifically for 
high contrast imaging, these instruments have the potential to 
probe dust-yatter ed light at small inner working angles (e.g., 
iPerrin et H] 120141) . With superior sensitivity and resolution, 
ALMA is capable for detailed mapping of the CO gas in this 
system. In addition, ALMA observations would trace the dis¬ 
tribution of larger grains which are more tightly coupled to the 
locations of the large parent bodies. Euture observations with 
these facilities hold great promise in further characterizing the 
dust distribution and dynamics. 


6. SUMMARY 

HD 131835 shows strong infrared excess, and here we 
present the discovery of the resolved debris disk in the 
mid-infrared. The debris disk’s properties can be constrained 
using all the available observations on HD 131835, including 
15 photometry points, the IRS spectrum, and resolved images 
at 11.7 /im and 18.3 /im. Erom our modeling result, the 
disk clearly can not be described by a single continuous 
population of modified blackbody grains. The images alone 
can be described by a two grain population model which is 
composed of a continuous power-law disk and a ring. The 
continuous component contains hotter grains and dominates 
the emission in the 11.7 jim image whereas the ring com¬ 
ponent contains cooler grains and dominates the emission 
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at 18.3 jim image. However, in order to obtain a good fit 
to the SED simultaneously, an additional ring component is 
needed. This three-component model is composed of a hot 
continuous power-law disk, a warm ring, and a cold ring. In 
this model, the disk fiuxes in the imaged wavelengths are 
contributed primarily from the hot continuous disk and the 
warm ring. Since the cold ring peaks at a longer wavelength, 
this third component does not show up in the mid-infrared 
images so its spatial location is relatively unconstrained. The 
excess emission in far infrared and submillimeter, however, 
can be well described by this third component. Starting 


with the simplest model, we found that a model with three 
components and therefore three grain populations can well 
describe the images and the SED simultaneously. 
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